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ABSTRACT 

Time dependent variable measures were obtained for all 
major aircraft accidents between July 1971 and July 197^. 
Using these time dependent variables and functional forms 
of these variables, a regression analysis was performed 
for each of eight major aviation commands. By using these 
functional forms of the variables, a relatively high amount 
of variance in aircraft accident rate was accounted for at 
a high confidence level in some commands. When reviewing 
the results of the eight major commands- considered, it was 
particularly noted that the variables most instrumental 
in explaining the variance in aircraft accident rate were 
not all pilot oriented but were variables interpreted as 
being related either to pilot experience level, pilot 
proficiency or aircraft condition. 
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I. INTRODUCTION 



Aviation safety is one of the major tools used to obtain 
a Navy goal of fleet readiness. The reports necessary to 
provide the Naval Safety Center with a means by which to 
judge the Navy's progress in aviation safety are delineated 
in OPNAVINST 3750*6 (Series). Through analysis of reported 
mishaps, the safety center has the ability to determine 
correctable material and personnel deficiencies. Through 
the administrative effort expended in submitting reports, 
the Navy can be repaid many times in the savings of lives 
and material. 

Due to extensive post aircraft accident analysis by 
accident investigation teams, aircraft accidents have been 
broadly categorized in terms of causal factors. Causal 
factors are not to be confused with environmental factors 
but are any event, act, failure/malfunction, circumstance 
or occurrence, the presence or absence of which caused the 
accident. Contributing to accidents are poorly designed 
equipment, supervisory error, improper aircrew training, 
failure of a facility to provide adequate runway drainage, 
maintenance induced equipment malfunctions, etc. 

An accident is designated as a major accident if: 

(1) loss of life occurs; (2) complete loss of an aircraft 
is involved; or (3) substantial damage occurs to any 
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aircraft involved where substantial damage is defined in 
Appendix A of OPNAVINST 3750.6 (Series). The most common 
cause cited in the literature on aircraft accidents is 
pilot error. In a study by Brictson, et al. /~1969_7» 
approximately seventy-eight percent of the accidents stud- 
ied attributed pilot error as being the primary causal 
factor, with only eight percent being attributed to 
support personnel errors, and the remaining thirteen per- 
cent being attributed to aircraft failure, equipment failure 
and/or weather. 

A study conducted for the Royal Air Force by Goorney 
Z~1965_7» concluded that fatigue, personal worries, com- 
placency, emotional stress and lack of current flying expe- 
rience, directly contributed to pilot error. Goorney' s 
study states that excessive ground duties prior to flying 
attributed to fatigue and that personal worries concerning 
marriages, dating, housing, financial and work problems 
attributed to emotional stress. 

Collicot, et al. /~1972__7 categorized accident causal 
factors into pilot error, material failure, maintenance 
error, and miscellaneous other causes. The authors compared 
Air Force F-4 accidents with Navy-Marine F-4 accidents and 
attributed maintenance error disparities to the fact that 
Air Force F-4 aircraft only realized one- tenth of the can- 
nibalization of parts that Navy -Marine aircraft realized. 

The authors also singled out that the Navy-Marine officer 
job rotational policies resulted in lower in-type flight 
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hours for Navy -Marine aviators and thus adversely affected 
their pilot proficiency. The authors also noted that, in 
comparing dual piloted to single piloted aircraft, the dual 
piloted aircraft had fewer accidents per ten thousand hours 
than the single piloted aircraft. The authors conclusion 
was that pilot mental overload was critical in the determina- 
tion of a prime factor in pilot error. 

As a result of the Collicot study and a study by Kowalsky, 
et al. /~197^_7» "the emphasis that pilot induced accidents 
were caused by lack of pilot proficiency was shifted to 
implications that pilot error accidents may be attributed 
to a state of temporary mental overload where the pilots 
incorrectly evaluated information inputed to them during the 
overload period. 

Since the advent of the Naval Safety Center, a great deal 
of effort has been expended in maintaining extensive data 
banks of accident related information. Statistical Analysis 
of available data should enable predictive mathematical 
models to be constructed, as much of the data available can 
be construed as measures of pilot proficiency. Myers /~197^_7 
attempted such a statistical analysis, hypothesizing that 
measures of pilot proficiency and experience would suffice 
to form an adequate foundation for accident rate analysis. 

He applied statistical techniques of principle component 
analysis and cluster analysis to ten variables obtained from 
the Individual Flight Activity Reporting System (IFARS) which 
are submitted to the Safety Center. Comparing two groups of 
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pilots — one group being involved in aircraft accidents and 
the other group being accident free — his results were not 
as pronounced as desired. The fact that each group con- 
tained only fifty subjects is reason to believe that this 
small sample size suppressed his accident predictive results. 

Maxwell and Stuck! /~1975_7 applied regression analysis 
techniques to an all Navy study of aircraft accidents for 
the period July 1968 to June 197^* Their analysis indicated 
that k-S.Sfo of the total accident rate variability is explained 
by pilot related variable measures if one wishes to accept 
a statistical confidence level of 75% as meaningful. Results 
of their study corroborated with the previous studies cited, 
that pilot error of one sort or another is the single largest 
cause of aircraft accidents. 

The author of this writing agrees with Maxwell and Stuck! 
that sufficient data should be currently available, from 
which predictive capability is extractable. The variable 
nature of the monthly accident rate suggests underlying 
factors, causal and thus definable in their role of accident 
perpetration. Using statistical analysis to isolate variable 
measures associated with pilot proficiency and aircraft 
maintenance which vary either directly or indirectly with 
accident rate, then predictive and thus preventive know- 
ledge can assist in lowering the loss of human life and 
the dollar loss resulting from aircraft accidents. 
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II . PROBLEM DEFINITION 



The Naval Safety Center (NSC) , Norfolk, Virginia, has 
provided the author with major accident rates computed for 
aircraft accidents for eight major commands: (1) Commander 

Naval Air Forces Atlantic (COMNAVAIRLANT) ; (2) Fleet Marine 
Forces Atlantic (MARLANT) ; (3) Commander Naval Air Forces 
Pacific (COMNAVAIRPAC) ; (4) Fleet Marine Forces Pacific 
(MARPAC) ; (5) Commander Naval Air Training Command (CNATRA) ; 
(6) Commander Marine Training Command (MA.RTC) ; (?) Naval 
Reserves; (8) Commander Research Development Test and 
Evaluation/Naval Air Systems Command (RDT&E/NASC) ; for 
fiscal years 19^9 through 197^- The accident rate is defined 
as the total number of accidents in a given month multiplied 
by a constant factor of ten thousand and then divided by 
the total monthly hours flown for a given command. Major 
accidents by definition, as delineated in OPNAVINST 3750.6 
(Series), are characterized by loss of life, loss of aircraft, 
or extensive aircraft damage measured in necessary man-hours 
to effect repair, 

A variety of approaches to look for consistent trends 
or cyclic phenomena in aircraft accident rate of Navy-Marine 
accidents was undertaken by Poock /~1976_7. None were found 
that could pass any statistical tests to verify their 
existence beyond reasonable elements of chance . 
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The purpose of this paper is to further explore accident 
rate dependence on time related variable measures for the 
eight major commands listed above in hopes that one or more 
of these measures can be identified for later use in the 
reduction of accident rates and thus the reduction in the 
loss of lives and money. 
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Ill . ANALYTICAL PROCEDURES 



A. DATA SOURCE 

OPNAVINST 3750-6 (Series) promulgates the requirements 
and procedures for the reporting of each aircraft accident 
or incident involving all Navy and Marine aircraft. This 
establishes control over accident data with the goal of 
increasing aviation safety. The reports are forwarded 
to the Naval Safety Center (NSC) for inclusion in their 
master data bank. Accident data currently available from 
NSC can provide approximately eighty variable measures for 
each accident. 

B. DATA SELECTION 

The NSC data bank provides a ready source of data for 
each accident or data set. However, as this is a continuation 
of the work of Maxwell and Stuck! , the 2110 computer data 
cards obtained from NSC for their study were also used in 
this study. 

The initial step in the conduct of the current accident 
rate analysis was to select appropriate variable measures or 
data points. A data point for an accident was considered 
to be any suitable variable measure associated with the 
accident, and a data set consisted of data points for a 
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specific accident. Specific measures were chosen in mutual 
discussions with NSC personnel. 

Selection of appropriate data points required that each 
point be time dependent. Data point time dependency and 
subsequent selection was based on the variable descriptions 
contained in the Manual of Code Classification for Navy 
Aircraft Accident, Incident and Ground Accident Reporting 
(Code Manual) promulgated by NSC. 

A sufficient number of data sets had to be incorporated 
into the analysis to facilitate viable statistical results , 
but the span of time defined by the data sets had to be 
chosen with care. Unfortunately, of the six fiscal years 
of data available to the author, only fiscal years 1972, 

1973 » and 197^ were used in this study. Fiscal years 1969, 
1970 , and 1971 were eliminated primarily because some data 
sets were either incomplete or entirely void of information. 

A change in reporting procedures during 1972 alleviated this 
problem. Also, in the fiscal years chosen for the analysis, 
a particular accident too often had to be disregarded due to 
incomplete data or to the fact that the particular occurrence 
involved another aircraft and the other aircraft was deemed 
at fault. Therefore, either due to insufficient data or 
to an occurrence and not an accident, a total of 636 of 
the original 2110 data sets available were used. Data for 
each variable of these 636 data sets were complete and there 
is absolutely no data missing in the data sets used in 
the analysis. 
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Ten data points from each set of the original thirty 
data points requested by Maxwell and Stuck! were selected 
for inclusion in the current analysis. Nine of these data 
points are the same data points used in the Maxwell and 
Stuck! study with the addition of the data point "Years 
Experience as a Designated Naval Aviation" and the data 
point "Flight Purpose Code" deleted. In addition, the data 
points selected also allowed the author to construct four 
additional variables deemed as measures of experience and 
pilot proficiency. Table I lists the basic variables used 
in the study. The first ten variables listed are from the 
Naval Safety Center and variables eleven through fourteen 
are the constructed variables. 

Pilot age, years experience as a designated Naval aviator, 
experience, age at designation as a Naval aviator and total 
flight time in the aircraft involved in the reported accident 
have been considered as measures of pilot experience. If 
these variables are true measures of experience they should 
exhibit negative correlation with accident rate, assuming 
a relationship exists. 

There is the age old adage of "practice makes perfect" 
which this author construes as a measure of pilot proficiency. 
Included in the category of pilot proficiency were total 
flight time during the preceding ninety days, total night 

flight time during the preceding ninety nights, total day 

flight time during the preceding ninety days, the member 

of carrier landings in the last thirty days, the number of 

carrier landings in the last thirty nights, and the ratio 
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TABLE I 



DATA SET INCLUDED IN CURRENT STUDY 

1. Accident rate by month (RATE) 

2. Pilot's age (AGE) 

3 . Years experience as a designated Naval aviator (DNA) 

k. Total flight time in accident involved aircraft model 
(TTIME) 

5 . Total flight time during preceding ninety days (T0T90) 

6. Total night flight time during preceding ninety nights 
(NITE90) 

7 . Daylight carrier landings during preceding thirty days 
(CLDAY) 

8. Night carrier landings during preceding thirty nights 
(CLNIT) 

9. Number of aircraft tours (ACTOUR) 

10. Aircraft flight hours since last major or minor 
inspection (ACHRS) 

11. The ratio of years experience as a designated Naval 
aviator aund total flight time in the accident involved 
aircraft model (EXPER=DNA/TTIME) 

12. Age at designation as a Naval aviator (WINGS=AGE-DNA) 

13 . Total flight time during preceding ninety days minus 
total night flight time during preceding ninety nights 
( DAY90=TOT90 -NITE 90 ) 

14. The ratio of total night flight time during the pre- 
ceding ninety nights to total day flight time during 
the preceding ninety days (NITEDAY=NITE90/DAY90) 



16 



of night flight time to day flight time during the preceding 
ninety days and nights. 

To measure airframe age and general condition, two var- 
iables were selected. All Navy/Marine aircraft are required 
to undergo a Periodic Aircraft Rework (PAR) cycle for analysis 
and repair after having accumulated a specific number of 
flight hours. Thus, the number of aircraft tours was chosen 
as one variable. As a measure of aircraft condition, the 
aircraft flight hours since the last major or minor inspection 
was also selected as a variable . 

C. PRELIMINARY DATA PREPARATION 

A common theoretical proposition in parametric statistics 
states that changes in one variable can be explained by 
reference to changes in several other variables. Such a 
relationship is described in a simple way by a multiple 
linear regression equation as described in Appendix C . To 
use this statistical method the relationships between the 
variables require the following assumptions: (1) normality; 

(2) data must be in the interval measurement scale; 

(3) the number of observations exceeds the number of coef- 
ficients to be estimated; and (4) no exact linear relation- 
ships exist between any of the explanatory variables. 

Normality was assumed by invoking the Central Limit Theorem ' 
and the technique of averaging by month was used to transform 
the data into interval data. 
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Raw data for each of the first ten variables of Table I 
was averaged by month for each of the thirty-six months within 
the time span selected and was also used to compute the values 
for the remaining variables of Table I. 

D. THE ANALYSIS TECHNIQUE 

To conduct the statistical analysis of the data set the • 
author selected the forward (stepwise) multiple regression 
computer program package developed by Jae-On Kim and Frank 
J. Kahout at the University of Iowa. The program is included 
in the Statistical Package for the Social Sciences (SPSS) 
compiled and edited by Nie, et al. /~1975__7. 

Multiple regression may be viewed as a descriptive tool 
whereby the linear dependence of one variable on other var- 
iables is summarized and decomposed. Kim and Kahout state 
that the most important uses of the technique as a descriptive 
tool are: 

(1) to find the best linear prediction equation 
and evaluate its prediction accuracy; (2) to 
control for other confounding factors in order 
to evaluate the contribution of a specific 
variable or set of variables; and (3) to find 
structural relations and provide explanations 
for seemingly complex multivariate relation- 
ships, such as is done in path analysis. 

However, the primary purpose is to evaluate and measure over- 
all dependence of a specified (dependent) variable on a set 
of other variables ( independent) . The dependent variable 
used for the current study was accident rate and the inde- 
pendent variables consisted of those variables listed as 2 
through l^’ in Table I and all the variables listed in Table II . 
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To explain the variables listed in Table II it must 
again be stated that regression analysis assumes the under- 
lying relationships among the variables are linear and addi- 
tive. There are many occasions for which such simple linear 
models are inadequate . Early efforts by the author dealt 
entirely with first order variables but examination of 
scatterplots of the residuals, which are conceived as meas- 
ures of the error component, indicated lack of linearity and 
the presence of curvilinearity . Therefore, functional terms 
were added to the variable list as suggested by Kim and Kahout. 
These additional variables are listed in Table II. 

The SPSS computer program is designed to provide the user 
with a considerable number of control options. The listwise 
deletion option was chosen for use by the author as it is 
the most conservative and accurate of the options. 

The forward (stepwise) multiple regression technique 
(Appendix C) is particularly useful in studies of the current 
type. It is appropriate to enter independent variables one 
by one on the basis of a pre-determined statistical criteria. 
This procedure is used when a researcher desires to isolate 
a subset of available predictor variables that will yield an 
optimal prediction equation with as few terms as possible. 

Draper and Smith 19 6 6_7 conclude the following in their 

discussion of various regression procedures: "We believe 

this (stepwise regression) to be the best of the variable 
selection procedures discussed and recommend its use." 

The percentage of the accident rate variability explained 
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TABLE II 



1. 


CONSTRUCTED 

AGE2 = (AGE)^ 


VARIABLES 

20. 


CLDAYI = 


1.0/CLDAY2 


2. 


TTIME2 = (TTIME)^ 


21. 


CLNITI = 


1.0/CLNIT2 


3. 


T0T902 = (T0T90)^ 


22. 


ACTOURI 


= 1.0/ACT0UR2 


4. 


WINGS 2 = (WINGS) ^ 


23. 


ACHRSI = 


1.0/ACHRS2 


5. 


DAY902 = (DAY90)^ 


2 ^-. 


DNAI = 1.0/DNA2 


6 . 


CLDAY2 = (CLDAY)^ 


25. 


EXPERI = 


1 . 0/EXPER2 


7. 


CLNIT2 = (CLNIT)^ 


26 . 


WINGSI = 


1.0/WINGS2 


8. 


ACHRS2 = (ACHRS)^ 


27. 


RTAGE = 


(AGE)^ 

(DNA)2 

= (TTIME)2 
= (T0T90)2 


9. 


EXPER2 = (EXPER)^ 


28. 


RTDNA = 


10. 


DNA2 = (DNA)^ 


29. 


RTTTIME 


11. 


NITE902 = (NITE90)^ 


30. 


RTT0T90 


12. 


ACT0UR2 = (ACTOUR)^ 


31. 


RTDAY90- 


= (DAY90)2 


13. 


NITEDAY2 = (NITEDAY) ^ 


32. 


RTCLDAY 


= (CLDAY)^ 
= (CLNIT)2 
= (ACHRS)2 


14. 


NITEDAYI = 1.0/NITEDAY2 


33. 


RTCLNIT 


15. 


AGEI = 1.0/AGE2 


3^. 


RTACHRS 


16 . 


TTIMEI = 1.0/TTIME2 


35. 


RTEXPER 


= (EXPER)2 


17. 


T0T90I = 1.0/T0T902 


36. 


RTWINGS 


= (WINGS) 2 
= (NITE90)2 
= (ACT0UR)2 


18. 


NITE90I = 1.0/NITE902 


37. 


RTNITE90 


19. 


DAY90I = 1.0/DAY902 


38. 


RTACTOUR 






39. 


RTNITEDAY = (NITEDAY)^ 
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by the time related independent variables entered is exactly 
the type of statistical output required to ascertain the 
causal factors responsible for that variability. 

Thirteen separate regression calculations were attempted 
for each of the eight commands considered. Each of the 
separate regression calculations were permutations of the 
variables listed in Table I and the variables listed in 
Table II. Regression calculations using only the first 
eight variables of Table I or regression calculations using 
the first ten variables' of Table I led to the suspicion of 
curvilinear ity in the data as a result of an inspection of the 
residual scatter plots. Therefore, the author attempted 
various regression calculations using combinations of the 
variables of Table I with either the squares of those var- 
iables, the inverse squares of those variables, the square 
roots of those variables or any combinations thereof. 

There is no unique procedure for selecting the best 
criteria and personal judgement is a necessary part of any 
decision reached by the analyst. The author set up guide- 
lines to choose the "best” equation as follows: (1) the 

regression equation would be limited to five or six inde- 
pendent variables; (2) the regression equation chosen should 
account for the highest amount of the variance in aircraft 
accident rate; (3) the residual plots would be examined 
for the best indication of minimum residuals; and (4) the 
statistical significance of each regression would be cal- 
culated as outlined in Appendix D and the regression equation 
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yielding the greatest amount of accountable variance with 
the highest confidence level would be the equation selected. 



22 



IV. RESULTS 



A. COMNAVAIRLANT 

Early efforts at regression analysis by the author used 
only the variables listed in Table I. These efforts led to 
the suspicion of curvilinearity as a result of examination of 
the residual plots. The order of variable inclusion at the 
early stage was: (1) total flight time during the ninety 

days preceding the accident (T0T90) ; and (2) aircraft flight 
hours since the last major or minor inspection (ACHRS) . In 
equation form the regression became: 

RATE = 0.35295 + 0 ,00658(T0T90) + 0 . 00l63( ACHRS) 

This equation accounted for 19-5?^ of the variance in aircraft 
accident rate at a 95?^ confidence level. 

Final regression results using functional forms of the 
variables indicated that the hierarchical order of variable 
inclusion as governed by the individual variable contributions 
towards explaining accident rate variance was: (1) the square 

root of aircraft flight hours since the last, major or minor 
inspection (RTACHRS) ; (2) the square of aircraft flight 
hours since the last major or minor inspection (ACHRS2) ; 

(3) the inverse square of the total flight time during the 
ninety days preceding the accident (TOT 90 I) j (4) the inverse 
square of the quotient years experience as a designated 
Haval aviator and total flight time in the accident involved 
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aircraft model (EXPERI) ; and (5) the square root of the 
quotient years experience as a designated Naval aviator 
and total flight time in the accident involved aircraft 
model (RTEXPER) . Table III lists the simple correlation 
coefficients and Table IV is a summary listing of computer 
output provided by the SPSS package. In equation form the 
regression became; 

RATE = 0.89771 + 0.118 23 (RTACHRS) - 0 .00002(ACHRS2) 

- 464.21653(T0T90I) - 0 .00002(EXPERI) 

- 4. 7699 6 (RTEXPER) 

or in terms of Table I variables 

RATE = 0.89771 + 0 .11823(ACHRS)^ - 0 . 00002( ACHRS) ^ 

- 464 . 21653 ( T0T90 ) - 0.00002(TTIME/DNA)^ 

- 4.76996(DNA/TTIME)2 

These equations account for 42.9% of the variance in aircraft 
accident rate at a 99% confidence level. Therefore, by 
using functional forms of the variables, an additional 
23 . 4 ^ of the variance was accounted for — a substantial 
improvement in the results. 

The net effect of the functional forms of hours since 
the last major or minor inspection is positive. The remain- 
ing variables are interpreted as being measures of either 
pilot experience or pilot proficiency and they have a net 
effect of reducing aircraft accident rate. 
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B. lyiAJRLANT 



Using the same early methods as those used for GOMNAVAIRLANT 
the order of variable inclusion at the early stage was: 

(1) years experience as a designated Naval aviator (DNA) ; 

(2) daylight carrier landings during the preceding thirty 

days (CLDAY) ; (3) "the quotient of years experience as a 
designated Naval aviator and total flight time in the accident 
involved aircraft model (EXPER) ; (4) the quotient of night 

flight hours to day flight hours during the preceding ninety 
days (NITEDAY) ; and (5) total flight time in accident 
involved aircraft model. In equation form the regression 
became: 

RATE = 0.96753 + 0.09865(DNA) + 0 .08248(CLDAY) 

- 20. 27987 (EXPER) + 1 . 41098 (NITEDAY) 

- 0 .00048 (TTIME) 

This equation accounted for 32 . 6 ^ of the variance in air- 
craft accident rate at a 75^ confidence level. 

Final regression results using functions of the var- 
iables indicated that the hierarchical order of variable 
inclusion as governed by the individual variable contri- 
butions towards explaining accident rate variance was : 

(1) the square root of years experience as a designated 
Naval aviator (RTDNA) ; (2) the square root of total 

night flight time during the preceding ninety nights 
(RTNITE 90 ); (3) total night flight time during the 
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preceding ninety nights (NITE 90 ) ; (4) the quotient of 

years experience as a designated Naval aviator and total 
flight time in the accident involved aircraft model (EXPER) ; 
( 5 ) total flight time in the accident involved aircraft 
model (TTIME ) 5 and ( 6 ) the square of night carrier land- 
ings during the preceding thirty nights (C1NIT2) . Table V 
lists the simple correlation coefficients and Table VI is a 
summary listing of computer output provided by the SPSS 
package. In equation form the regression became: 

RATE = - 0.17128 + 0.88 657 (RTDNA) + 0 . 52634(RTNITE90) 

-0.06351(NITE90) - 57 . 3 70 27 (EXPER) - 0. 00 115( TTIME) 
+ 0.21877(CLNIT2) 

or in terms of Table I variables 

i i 

RATE = - 0.17128 + 0.88657(DNA)2 + 0 . 52634(NITE90) ® 

- 0 . 0635 KNITE 90 ) - 57 . 370 27 (DNA/TTIME) 

- 0.00115( TTIME) + 0.21877(CLNIT)^ 

These equations account for 5^.5?^ of "the variance in air- 
craft accident rate at a 99 ?^ confidence level versus the 
32 . 6 ^ accounted for initially. A substantial improvement 
was again observed when using functional forms of the 
variables in the regression analysis. 

When reviewing the correlation matrix (Table V) it is 
noted that every variable present in the equation is posi- 
tively correlated with accident rate except for EXPER (a 
measure of pilot experience) and the square of night 
carrier ’ landings in the last thirty days (a measure of pilot 
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proficiency) . The net affect of the functional forms of 
night flight hours during the last ninety nights is positive, 
an indication that night flying may he a hazardous evolution 
and add to accident rate. 
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MARLANT MATRIX OF SIMPLE 
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(CONSTANT) -0.17128 



C. COMNAVAIRPAC 



Initial efforts yielded the following order of variable 
inclusion: (1) total flight time during the preceding 

ninety days (DAY90) ; (2) the number of aircraft tours 

(ACTOUR) ; and (3) years experience as a designated Naval 
aviator. In equation form the regression became: 

RATE = -0.079^^+ 0.01289(DAY90) + 0 . 1056l(ACT0UR) 

+ 0.04137(DNA) 

This equation accounted for 33*8?^ of the variance in air- 
craft accident rate at a 99 % confidence level. 

Final regression results using functions of the var- 
iables indicated that the hierarchical order of variable 
inclusion as governed by the individual variable contri- 
butions towards explaining accident rate variance was : 

(1) the inverse square of total day flight time during 
the preceding ninety days (DAY90I); (2) the square of 

the age of the pilot at designation as a Naval aviator 
(WINGS2) ; (3) the number of aircraft tours (ACTOUR); 

(4) the inverse square of aircraft hours since the last 
major or minor inspection (ACHRSI); and (5) the square 
of the hours since the last major or minor inspection 
(ACHRS2) . Table VII lists the simple correlation coef- 
ficients and Table VIII is a summary listing of computer 
output provided by the SPSS package. In equation form the 
regression became: 
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RATE = 1.47479 - 670.77272(DAY90I) - 0 .00134(WINGS2) 
+ 0. 15733 (ACTOUR) + 304.94442(ACHRSI) 



+ 0.00001(ACHRS2) 

or in terras of Table I variables 

RATE = 1.47479 - 670.77272/(DAY90)^ - 0.00134(AGE - DNA)^ 

+ 0.15733(ACT0UR) + 304.94442(ACHRS)^ 

+ O.OOOOl(ACHRS)^ 

These equations account for 45. 7?^ of the variance in air- 
craft accident rate at a 99-5^ confidence level, an increase 
of 11 . 9 ?^ in accounting for variance. 

It is noted in these equations that the nuraber of 
aircraft tours and functional forras of hours since the 
last raajor or minor inspection have a positive effect on 
accident rate. Conversely, variables interpreted as 
measures of pilot proficiency and experience level have 
a negative effect. 
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COMNAVAIRPAC MATRIX OF 
SIMPLE CORRELATION COEFFICIENTS 
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D . MARPAC 



Initial efforts at regression analysis yielded very poor 
results for this command. The only variable to be included 
in the equation was the quotient of night flight hours to 
day flight hours during the preceding ninety days . In 
equation form the regression became: 

RATE = 1.75337 - 2.47720(NITE90/DAY90) 

This equation accounted for 9 • of the variance in air- 
craft accident rate at a 75?^ confidence level. 

Further regression efforts using functions of the var- 
iables led to much better results with the hierarchical 
order of variable inclusion being: (1) the square of 

the number of aircraft tours (ACT0UR2) ; (2) the number 

of aircraft tours (ACTOUR) : (3) "the square of total 

night flight time during the preceding ninety nights 
(NITE902) ; (4) the square of the quotient total night 

flight hours to total day flight hours during the preceding 
ninety days (NITEDAY2) ; (5) aircraft flight hours since 

the last major or minor inspection (ACHRS) ; and (6) the 
square of age at designation as a Naval aviator (WINGS2) . 
Table IX lists the simple correlation coefficients and 
Table X is a summary listing of computer output provided 
by the SPSS package. In equation form the regression 
became : 
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RATE = 1.26335 - 0.32047(ACT0UR2) + 2. 21178(ACT0UR) 

+ 0 .01042(NITE902) - 19.20769(NITEDAY2) 

- 0. 00565 (ACHRS) - 0.00^62(WINGS2) 

or in terms of Table I variables 

RATE = 1.26335 - 0.32047(ACT0UR)^ + 2 . 21178( ACTOUR) 

+ 0.01042(NITE90)^ - 19 . 20769 (NITE90/DAY90) ^ 

- 0.00565(ACHRS) - 0.00462(AGE-DNA)^ 

Either equation accounts for 39 • of the variance in air- 
craft accident rate at a 95 % confidence level, a substantial 
improvement over the 9 <2% of variance accounted for initially. 

Reviewing the correlation matrix it is noted that every 
variable present in the equation is negatively correlated 
with accident rate. However, the intercorrelations of the 
variables has the affect of causing some of the signs of 
the coefficients in the equation to change, with the net 
effect of the functional forms of aircraft tours being 
positive. 
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MARPAC MATRIX OF SIMPLE 
CORRELATION COEFFICIENTS 
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ACHRS 0.6^657 0.41805 0.32852 - 0.00565 - 0.39225 
WINGS2 0.70225 0.49315 0.39178 -0.00462 - 0.30796 
(CONSTANT) 1.26335 



E. CNATRA 



All attempts at regression analysis for this command 
yielded essentially the same results. Early efforts using 
only the variables listed in Table I yielded the following 
order of variable inclusion: (1) the quotient of total 

night flight hours to total day flight hours during the 
preceding ninety days (NITEDAY) ; (2) total flight time 

in accident involved aircraft model (TTIME) ; (3) the 

quotient of total flight time in accident involved aircraft 
model to years experience as a designated Naval aviator 
(TTIME/DNA) ; (4) aircraft flight hours since the last 

major or minor inspection (ACHRS ) ; and (5) the quotient 
of years experience as a designated Naval aviator to total 
flight time in accident involved aircraft model. In equation 
form the regression became: 

RATE = 5.22976 + 51.98496(NITE90/DAY90) + 0 .00955( TTIME) 

- 0.01973(TTIME/DNA) - 0 .04198 (ACHRS) 

- 216.41848(DNA/TTIME) 

This equation accounted for 32. 3?^ of the variance in air- 
craft accident rate at a 90 ^ confidence level. 

Later regression results using functions of the variables 
yielded: (1) the square root of aircraft flight hours since 

the last major or minor accident (RTACHRS) ; (2) the square 

of aircraft hours since the last major or minor inspection 
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(ACHRS) j (3) 'the square of the total day flight hours 
flown during the preceding ninety days (DAY902) ; (4) the 

square root of the quotient of total flight time in the 
accident involved aircraft model to years experience as a 
designated Naval aviator ( 1 . O/RTEXPER) ; and (5) the 
quotient of total flight time in accident involved air- 
craft model and years experience as a designated Naval 
aviator. Table XI lists the simple correlation coefficients 
and Table XII is a summary listing of computer output pro- 
vided by the SPSS package. In equation form the regression 
became : 

RATE = -0.31346+ 0.11 248 (RTACHRS) - 0 . 00003( ACHRS2) 

- 0.00002(DAY902) + 0 . 03482( 1. O/RTEXPER) 

- 0.00092(1.0/EXPER) 

or in terms of Table I variables 

RATE = - 0.31346 + 0 .11248 (ACHRS - 0 . 00003(ACHRS) ^ 

O 1 . 

- 0 .00002(DAY90) + 0 . 03482 (TTIME/DNA ) ^ 

- 0 .00092(TTIME/DNA) 

Either equation accounts for 33*0^ of the variance in air- 
craft rate at a confidence level of 90%; representing no 
significant change in the results. The functional forms 
of aircraft hours since the last major or minor inspection 
and the reciprocal foimis of EXPER have a net positive 
effect while the square of day flight time during the last 
ninety days has a negative affect. 
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CNATRA MATRIX OF SIMPLE 
CORRELATION COEFFICIENTS 
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F . NAVAL RESERVES 



Initial efforts yielded the following order of variable 
inclusion: ( 1 ) total night flight time during the preceding 

ninety nights (NITE 90 ); ( 2 ) the quotient of total night 
flight time and total day flight time during the preceding 
ninety days and nights (NITE90/DAY90) ; (3) total flight 

time during the ninety days preceding the accident (T0T90); 
and (4) the number of daylight carrier landings during the 
preceding thirty days. In equation form the regression 
became : 

RATE = 0.97395 + 0.06986(NITE90) - 1. 2295^(NITE90/DAY90) 

- 0.00777(T0T90) - 0. 033 ^ 3 (CLDAY) 

This equation accounted for 71*^ of the variance in air- 
craft accident rate at a 99^ confidence level. 

Final regression results using functions of the variables 
indicated that the hierarchical order of variable inclusion 
as governed by the individual variable contributions towards 
explaining aircraft accident rate variance was: (1) the 

square of total night flight time during the preceding 
ninety nights (NITE902) j (2) the square of the quotient 
of total night flight time and total day flight time during 
the preceding ninety days (NITEDAY2) ; (3) "the inverse 

square of aircraft flight hours since the last major or 
minor inspection (ACHRSI) ; (4) the square root of daylight 

carrier landings during the preceding thirty days: and 
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(5) "the inverse square of the quotient of years experience 
as a designated Naval aviator and total flight time in the 
accident involved aircraft model (EXPERI) . Table XIII lists 
the simple correlation coefficients and Table XIV is a summary 
listing of computer output provided by the SPSS package. In 
equation formi the regression became: 

RATE = 0.83553 + 0 .00084(NITE902) - 0 . 56367(NITEDAY2) 

- 311 . 23233 (ACHRSI) - 0. 22889 (RTCLDAY) 

+ 0. 0000 2( EXPERI) 

or in terms of Table I variables 

RATE = 0.83553 + 0.00084(NITE90)^ - 0 . 56367 (NITE90/DAY90) ^ 

- 311- 23233 (l.O/ACHRS)^ - 0.22889(CLDAY)2 

+ 0.00002(TTIME/DNA)^ 

Either equation accounts for 81. 3?^ of the variance in air- 
craft accident rate at a 99 •9?^ confidence level, a 9 . 9 % 
increase of accountable variance over initial efforts. The 
net effect of the square of night flight hours during the 
preceding ninety nights and the quotient of night flight 
hours and day flight hours during the preceding ninety 
nights and days is negative. This may indicate that the 
higher the ratio of night flight hours to day flight hours 
(a measure of pilot proficiency), the lower the accident 
rate . 
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NAVAL RESERVES MATRIX OF 
SIMPLE CORRELATION COEFFICIENTS 
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G . MARTC 



Initial efforts using only the variables of Table I 
yielded the following order of variable inclusion: (1) the 

quotient of total night flight hours and total day flight 
hours during the preceding ninety days (NITEDAY) ; (2) total 

flight time in accident involved aircraft model (TTIME) ; 

(3) the quotient of total flight time in accident involved 
aircraft model and years designated as a Naval aviator 
(TTIME/DNA) ; (4) age at designation as a Naval aviator 

(WINGS ) f and (5) number of aircraft tours (ACTOUR) . In 
equation form the regression became: 

rate = - 45.84436 + 32 . 239 27 (NITEDAY) + 0 . 0270? (TTIME) 

- 0.15004(1.0/EXPER) + 2. 20096 (WINGS) 

- 0 .75295(ACT0UR) 

This equation accounted for 68.0fo of the variance in air- 
craft accident rate at a 95^ confidence level. 

Final regression results using functions of the variables 
indicated that the hierarchical order of variable inclusion 
as governed by the individual variable contributions towards 
explaining aircraft accident rate variance was: (1) the 

square of the quotient total night flight time and total 
day flight time during the preceding ninety days (NITEDAY2) ; 
(2) total flight time in accident involved aircraft model 
(TTIME) ; (3) "the inverse square of the quotient years expe- 

rience as a designated Naval aviator and total flight time in 
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% 



I 






I 




the accident involved aircraft model ( EXPERT ) ; (4) total 

night flight time during the preceding ninety nights (NITE 90 ) ; 
3 Tid (5) the square of age at designation as a Naval aviator 
(WINGS2) . Table XV lists the simple correlation coefficients 
and Table XYI is a summary listing of computer output provided 
by the SPSS package. In equation form the regression became: 

RATE = - 0.36132 + 65 . 28325(NITEDAY2) + 0 . 00328(TTIME) 

- 0.00008(EXPERI) - 0 . 4l913(NITE90 ) 

+ 0.00385(WINGS2) 

or in terms of Table I variables 

RATE = - 0.36132 + 65 . 28325(NITE90/DAY90)^ + 0 .00328(TTIME) 

- 0.00008(TTIME/DNA)^ - 0 .41913(NITE90) 

+ 0.00385(AGE-DNA)^ 

Either equation accounts for 79.9% of the variance in air- 
craft accident rate at a 95 % confidence level, a 10 . 9 % increase 
of accountable variance over initial efforts. 

It should be noted for this command that only one variable 
interpreted as a measure of experience (EXPERT) is negatively 
correlated with accident rate. Another measure of experience 
(WINGS2) and three variables interpreted as measures of pilot 
proficiency are all positively correlated with accident rate. 
This positive correlation of pilot proficiency aind experience 
variables is contrary to what the author hypothesized. The 
study by Poock /~1976_7 also concluded that this command 
produced- unusual results. 



43 



MARTC MATRIX OF SIMPLE 
CORRELATION COEFFICIENTS 
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H. RDT&E/NASC 



These two commands were grouped together due to the fact 
that neither command had enough data points to be analyzed 
if considered as individual commands. Initial efforts using 
only the variables of Table I yielded the following order of 
variable inclusion: (1) total night flight hours during the 

preceding ninety nights (NITE90) ; (2) the quotient of 

years experience as a designated Naval aviator and total 
flight time in accident involved aircraft model (EXPER) ; 
and (3) pilot's age. In equation form the regression 
became : 

RATE = 0.93622 - 0 .08430 (NITE90) + 0 . 3l653( EXPER) 

- 0.01623(AGE) 

This equation accounted for > 5 % of the variance in air- 
craft accident rate at a 95^ confidence level. 

Final regression results using functions of the variables 
indicated that the hierarchical order of variable inclusion 
as governed by the individual variable contributions towards 
explaining aircraft accident rate variance was: (1) the 

square of total night flight during the preceding ninety 
nights (NITE 902) 5 (2) the square root of total flight 

time in accident involved aircraft model (RTTTIME) ; (3) the 

inverse square of years experience as a designated Naval 
aviator (DNAI) ; (4) the square root of the inverse quotient 
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of years experience as a designated Naval aviator and 
total flight time in the accident involved aircraft model 
(I.O/RTEXPER) } (5) the inverse square of total flight 

time during the preceding ninety days (T0T90I) ; and (6) the 
square root of total flight time during the preceding ninety 
days (RTT0T90) . Table XVII lists the simple correlation 
coefficients and Table XVIII is a summary listing of com- 
puter output provided by the SPSS package. In equation 
form the regression became; 

RATE = - 1.21331 + 0.71262(NITE902) + 0 . 23598(RTTTIME) 

+ I9.6852KDNAI) - 0. 61421 (l.O/RTEXPER) 

+ 182.91525(T0T90I) + 0 . 38527(RTT0T90) 

or in terms of Table I variables 

p i 

RATE = - 1.21331 + 0.71262(NITE90) + 0 . 23598(TTIME) ^ 

1 . 

+ I9.6852KI.O/DNA) - 0.61421(TTIME/DNA)2 
+ 182. 91525(1. 0/TOT90)^ + 0.38527(T0T90)^ 

Either equation accounts for 80.8?^ of the variance in air- 
craft accident rate at a 95?^ confidence level, a 25. 3?^ 
increase over initial efforts. 
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rdt&e/nasc matrix of simple 

CORRELATION COEFFICIENTS 
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V. DISCUSSION 



Reviewing the results it can be seen that as functional 
forms of the variables were introduced into the regression 
analysis that more variance in aircraft accident rate was 
accounted for at a higher confidence level. This was true 
for all eight commands considered. 

Observing the signs of the correlation coefficients 
listed in Table XIX, it can be seen that the variable EXPER, 
which is the ratio of years experience as a designated 
Naval aviator and total flight time in accident involved 
aircraft model, was negatively correlated with accident 
rate in all cases. This indicates that the more flight 
time in a particular aircraft model that a pilot possesses 
per year of designated Naval aviator service that he will 
become more safe and thus have fewer accidents. This could 
suggest that the more experienced a pilot becomes the larger 
a repertoire of near tragedies he has to draw from and thus 
the more reminders or analogies he has to compare with 
current situations. 

Conversely, the variable CLNIT, which is the number of 
night carrier landings a pilot has during the last thirty 
nights, showed positive correlation with aircraft accident 
rate in the three commands where the variable was present 
(COIVINAVAIRLANT, MARLANT, COMAVAIRPAC ) . This may indicate 
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SIGNS OF CORRELATION COEFFICIENTS 
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that night carrier landings are more taxing evolutions 
and tend to add to aircraft accident rate. 

It has been common belief that pilot proficiency lessens 
aircraft accident rate. However, when again reviewing the 
correlation coefficients, it is observed that the correlation 
of total flight time during the preceding ninety days (T0T90) 
and day flight time during the preceding ninety days (DAY90) 
is positively correlated with aircraft accident rate except 
in the cases of CNATRA and MARTC (both training commands). 

It is also observed that the variable night flight time 
during the preceding ninety nights is positively correlated 
with aircraft accident rate in all commands except CNATRA. 
Considering that the training commands are on a schedule 
where the pilot should not experience fatigue, the positive 
correlation of these variables with aircraft accident rate 
would tend to indicate that if flight time is not uniformly 
spread out that practice does not necessarily maJce perfect 
and that a possible fatigue factor may be causing an increase 
in aircraft accident rate. Fatigue could be construed as 
having the majority of flight time during the preceding 
ninety days massed into a short period of time instead of 
being distributed uniformly over the ninety day period. (It 
is understood that data for a thirty day period is now being 
recorded.) Another explanation would support Goorney's 
supposition that pilot complacency may increase directly 
as the number of hours flown and thus contribute to aircraft 
accidents. In this case, new pilots in the training command 
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most likely will not become complacent or they will not 
receive their wings. 

The remaining variables listed in Table XIX did not 
indicate a dominant sign of the correlation coefficient 
common to all commands. Therefore, the author will not 
attempt to explain the remaining variables individually. 

When reviewing the results of the major commands con- 
sidered, it should be noted that the variables most instru- 
mental in explaining the variance in aircraft accident 
rate are not all pilot oriented variables as was shown in 
the Maxwell and Stuck! study which was an overall Navy/ 

Marine study. The overall results don’t seem to apply 
when taking a more microscopic look at individual commands. 
COMNAVAIRLANT, COMNAVAIRPAC , MARPAC, GNATRA and Naval Reserves 
regression equations all contain variables interpreted as 
being either related to experience level, pilot proficiency, 
or aircraft condition. Only MARLANT, MARTC, and RDT&E/NASC 
accoiint for the variance in aircraft accident rate by 
variables interpreted as being pilot oriented variables 
in their respective regression equations. Condition of 
the aircraft must therefore also be considered when attempt- 
ing to account for the variance in aircraft accident rate. 
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VI. RECOMMENDATIONS 



The current study includes eight major commands in its 
treatment of aircraft accident rate. Functional forms of 
the variables produced the best results. For this type of 
study however, the author recommends that the distributions 
for each variable be ascertained and that these distributions 
whether they be gamma, beta, exponential, etc. be used in 
an attempt to analyze aircraft accident rate. 

This study was limited to consideration of only accident 
involved pilots and accident involved aircraft. A future 
study could be done comparing these pilots and aircraft 
to accident free pilots and aircraft in order to ascertain^ 
whether these groups are from the same distribution. It 
would then be possible to ascertain critical points in a 
pilot's career or the life of an aircraft. 

A more microscopic study than this study may seem appro- 
priate. However, subdividing the data any further into 
Airwings may limit the analyst to only a small number of 
data points . 
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AVERAGE MONTHLY DATA POINTS 
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?4 o 2 0.73 41.00 17.50 



TABLE A4 (CONTINUED) 
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TABLE A 5 

AVERAGE MONTHLY DATA POINTS 
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NOTE: Accidents of July 1971 , January 1973 » and January 1974 were deleted due to 

incomplete data. 
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NAVAL RESERVES 



m 

PC, 

X 

o 

c 



oooooooooooooooooo cn 

00000000^0^0000000 CO 



0^r-lC^-r-^\r^O^OCM >4^00 o\o\^oono-^ on 

0 4j“ CO (TNCnvO H C\2 ^rH O £>- C}\ 

r— 1 CM r— 1 r— 1 CM r— 1 CM rH 



PC 

:d 

o 

E-» 

O 

c 



oooooooooooooooooo 

oooooooooooooooooo 



MDMD O-CMmD O-CM Cn O- 00 MD O- 



00 

CM 



oooooooooooooooooo 

oooooooooooooooooo 



o 

o 



o 



oooooooooooooooooo 



>-i 

c 

Q 

o 



oooooooooooooooooo C5\ 

oooooooooooooooooo CO 



OONOOOOOOOOOOOOOOOO o 



o 

CT\ 



w 






oooooooooooooooooo 

0000000M-^\r^00000M-^000 



4itO OOOM*^C>-cOiHCO OMD OC^O-^iH O-MD 
rH CM 



ON 



CO 



oooooooooooooooooo c>- 

OOOOOOO^OOM^OOO^OOO CTn 

CMCN-CMCMOONOMD OM^OO^j- C^vO f— I O CM 
cn CM NO (T\ CM rH ON4:t (TN C^CO CM iH CM 



OOOOOOOOOOOOOOOOOO o- 
OOOOOOOOOOM^OOOOOOO 4:t 



CM vn o NO O O 4:^ CM O- CM XO O- O O 
O O CJN4:t NO CM 00 O O iH C^NOO CO o rH 

(ON CM rH iHjj- O-M^rH O rH xr\ rH O- 

rH rH I — 1 rH r — I rH 



OOOOOOOOOOOOOOM-^oOO tN- 
OOOOOOOO^OM-^oOOt^-OOM-^ o 

CM CM cm O CTs rH O CM O-CJN^it^it CM O- O 

rHrH rH rHCMrHiH rHrHrHrH 



OOOOOOOOOOOOOOOOOO CM 
Wl OOOOOOOOOOM^OOOM-^oOO O- 

O 

cl NO O-NO 00 CM (ON C:NM-^o^\r^J:t CM O O-OO O O cn 
(T^cnCMCMC^CMCMC^ Cnjzh 



C 

§ 



E-» 

E-» 



O 

ON 

O 

E-t 



C 

PC 



CnC^NOCMrHrHONOOOO M^NO rH CM MO m rH ^ 
xo vn CN-MO OONONOCMC^xr^CMxr^ X4^NO M^NO 0-0- 00 



OOOOOOOrHrHOrHOOOCMOOO O 



O 

PC 



00 OrH4:tM^O-<ONO CM 0-00 C7\ CM (T%c^ MO 

OrHrHOOOOrHOOOOOOrHOOO 
rHrHrHCMCMCMCMCMC^C0C^C0C0(?0 ^ 

{>-0-0-0-0-0-0-0-0-0.0.£N.O-0.(N.O.O.O. 



O 

> 

C 



61 



NOTE: The accident in August 1972 was deleted due to incomplete data. 
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AVERAGE MONTHLY DATA POINTS 
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AIRCRAFT ACCIDENTS BY MONTH AND BY COMMAND 
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APPENDIX B (CONTINUED) 
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APPENDIX C 



FORWARD (STEPWISE) MULTIPLE REGRESSION 

The basic principles of regression analysis may be 
extended to situations involving two or more independent 
variables. The general mathematical form of the (unstand- 
ardized) regression is 



Y’ = A + + 3^X2 + . . . + 

where Y’ represents the estimated value of the dependent 
variable, A is the Y intercept, the are regression 
coefficients, and the X^ are the independent variables. 

It is assumed that this is a complete set of variables 
from which the equation is to be chosen and includes any 
functions such as squares and inverse squares thought to 
be desirable and necessary. 

The A and B^ coefficients are selected in such a way as 

to minimize the sum of squared residuals, ^(Y-Y')^« By 

minimizing the squared residuals, the regression technique 

maximizes the correlation between the actual dependent 

variable (Y) and the estimated dependent variable (Y') 

while the correlation between the independent variables 

and the residual values (Y-Y') are reduced. 

The proportion of variamce of Y explained, i.e., the 

goodness of fit of the regression equation can be evaluated 

by examining the square of the multiple correlation (R ) , 

2 

where R is calculated by: 
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res 




SS 



SS 



reg 





(Y - Y-) 



or 






variation in Y explained by the combined 
linear influence of the independent variables 



total variation in Y 



where SS is the total variation or sum of squares in Y, 

is the sum of squared residuals, and SS is the 
r0s ^ reg 

regression of squares. 

The forward (stepwise) selection procedure available in 
the Statistical Package for the Social Sciences (SPSS ver- 
sion 6) used for this study inserts variables in turn until 
the regression equation is satisfactory. Independent var- 
iables are entered only if they meet certain statistical 
criteria. The order of inclusion is determined by the 
respective contribution of each variable that explains the 
greatest amount of variance previously unexplained by th« 
variables already in the equation. The variable that does 
this is the next variable to be entered. 

This order of insertion is determined by using the 
partial correlation coefficient as a measure of the impor- 
tance of variables not yet in the equation. The, basic 
procedure is as follows. First select the X most correlated 
with Y (suppose it is and find the linear regression 

equation Y = f(X^) . Next find the partial correlation 
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coefficient of X.(jVl) Y (after allowance for X,), 

i.e., find the correlation between the residuals from the 

A 

regression Y = f(X^) and the residuals from another regres- 
sion X.=f.(X,). The X. with the highest partial correlation 
J J -J- J 

with Y is now selected and the process continues. As each 
variable is entered into the regression, the following 

values are examined; 

2 

1. R , the multiple correlation coefficient; 

2. The partial F-test value for the variable most 
recently entered, which shows whether the variable 
has taken up a significant amount of variation 
over that removed by variables previously entered 
in the regression. 

As soon as the partial F value related to the most recently 
entered variable becomes nonsignificant the process is 
terminated . 
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APPENDIX D 



STATISTICAL TESTS FOR SIGNIFICANCE 



Regression procedures per se may be categorized as de- 
scriptive statistics. Regression analysis is commonly per- 
formed on sample data from which the researcher is either 
interested in estimating population parameters from sample 
regression statistics or to testing statistical hypothesis 
about the population parameters and determining confidence 
limits for estimates in testing the hypothesis. 

The overall test- for goodness of fit of the regression 
equation uses statistical inference procedures to test the 
null hypothesis that the sample of observations being analyzed 
has been drawn from a population in which the multiple cor- 
relation is equal to zero. An equivalent way of stating 
the null hypothesis, Ho, is that the next variable to be 
added in a forward regression would not add significantly 
to the explained variance in the dependent variable, Y, 
already accounted for by variables included in the regression 
equation. The alternate hypothesis, H^, directly contradicts 
the null hypothesis. 

The test statistic employed for the overall test is 





F = 



(1-R^)/(N-K-1) 
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where SS is the sum of squares explained by the entire 

1/6 g 

regression equation, SS is the residual (unexplained) 

X c o 

sum of squares, K is the number of independent variables 

in the equation and N is the sample size. The F ratio is 

distributed approximately as the F distribution with degrees 

of freedom K and (N-K-1) . 

2 2 

Adjusted R is an' R statistic adjusted for the num- 
ber of independent variables in the equation and the number 
of cases. It is a more conservative estimate of the percent 

of variance explained, especially when the sample size is 

2 

small. The adjusted R formula uses unbiased estimates 
of the error variance and the total variance of Y in the 
population. The formula used by SPSS is: 

Adjusted R^ = R^ - — (1-R^) 

2 

Becaused Adjusted R is a conservative estimate of the per- 
cent of variance explained, the author used the following 
F statistic to determine confidence levels for the regression 
equations obtained for each command considered: 

(Adjusted R^)/K 

p = 

(1-Adjusted R^)/(N-K-1) 
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